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Abstract. The 0(3) sigma model in two spatial dimensions admits topological (Bo- 
gomol'nyi) lower bound on its energy. This paper proposes a lattice version of this system 
which maintains the Bogomol'nyi bound and allows the explicit construction of static soli- 
tons on the lattice. Numerical simulations show that these lattice solitons are unstable 
under small perturbations; in fact, their size changes linearly with time. 
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1. Introduction. 

The nonlinear 0(3) sigma model in (2+1) dimensions is a popular model in theoretical 
physics; the static system is integrable and of Bogomol'nyi type (all minimal energy 
solutions can be obtained by solving the Bogomol'nyi equations). As a result, one can 
explicitly write down soliton solutions of arbitrary degree in term of rational functions Pj; 
but the model is scale invariant and therefore, its solitons have no fixed size and so their 
stability is a central question. Under small perturbations they shrink towards infinitely 
tall spikes of zero width or may spread out, with this expansion continuing indefinitely. 
That this indeed happens is confirmed by numerical experments [0, |3|. General time- 
dependent solutions cannot be constructed explicitly, and so it is natural to investigate 
numerical evolution techniques which discretize the partial differential equations. 

Given a continuum field theory, there are many different lattice systems which re- 
duce to it in the continuum limit. In systems where there are topological configurations 
(instantons, monopoles, etc) one often has a Bogomol'nyi bound which is related to the 
stability of the topological objects in question. If the bound is maintained on the lattice, 
the topological objects will be well-behaved even when their size is comparable to the 
lattice spacing. Lattice versions of these systems are important for purposes of numerical 
computations but they have, generally, not preserved the Bogomol'nyi bound. Few years 
ago, Leese || discretized the (unmodified) 0(3) sigma model in (2 + 1) dimensions. He 
imposed radial symmetry, made the radial coordinate r discrete and found a reduced 
lattice system with Bogomol'nyi bound. But, although the topological lower bound can 
be attained, the minimum-energy configurations are not explicit. On the other hand, 
Ward H described a lattice version of this model with Bogomol'nyi bound, without any 
symmetry constraint. In this general case, however, the lower bound cannot be attained. 

The purpose of this paper is to present a lattice version of the 0(3) sigma model in 
two space dimensions, in which the Bogomol'nyi bound is maintained. The primary aim 
is not to simulate the continuum system, but rather to define an alternative lattice system 
with similar properties but more convenient to study numerically. Following Leese, only 
field configurations for which the energy density (and not necessarily the fields) is radially 
symmetric will be considered here, so that in effect one obtains a one-dimensional system 
and therefore, the construction of the discrete Bogomol'nyi equations is less complicated. 
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Solutions of these equations (which were obtained analytically) are then used as the basis 
for a numerical study of soliton stability. In fact we study the shrinking of lattice solitons 
and try to estimate analytically the rate of this shrinking for a soliton and 2-soliton 
configuration. It was shown |J that a collision of two solitons can induce their shrinking, 
therefore we can then follow one of these solitons and study its behaviour. However, due 
to the soliton movement it is difficult to be very quantitative about the soliton shrinking. 
For this reason it is convenient the restriction to axial symmetry. It is worth remarking 
that if radial symmetry is imposed on the fields |7j then the model is integrable, and an 
inverse scattering transformation exist. But this restriction requires that the obtained 
solutions have topological charge zero. 

The rest of this paper is arranged as follows. In the next section we describe initially, 
the familiar continuum 0(3) sigma model in (2 + 1) dimensions then reparametrize the 
fields in order to impose radial symmetry and, finally, discretize the model. In section 3 
we study the dynamics of the 2-soliton configuration at low shrinking velocities using the 
slow-motion approximation, make approximate analytic predictions of its behaviour, and 
compare these with numerical results, while in section 4 we investigate the properties of 
the lattice 0(3) solitons numerically. 

2. The Lattice 0(3) Sigma Model. 

Let us begin with a brief review of the continuum 0(3) sigma model in two space 
dimensions. The field is a unit 3-vector field on R 2 (i.e. a smooth function from R 2 
to the target space S 2 ), with the boundary condition <fr — ► <p as r — > oo in R 2 . Here 
(fro is some fixed point on the image sphere S 2 . Hence there are distinct topological 
sectors classified by an integer k (topological charge), which represents the number of 
times R 2 is wrapped around S 2 . Roughly speaking, k is the number of solitons. The 
potential energy of the field is E p = (87r) _1 J[(d x (fr) 2 + (d y (fr) 2 } dxdy and the appropriate 
Bogomol'nyi argument gives the bound E p > \k\. There are fields which attain this lower 
bound (such minimum-energy fields will be called solitons in what follows). Since E p 
is invariant under the scaling transformation (fr(x^) i— > <p(Xx^) these configurations are 
metastable rather than stable (their size is not fixed). 

From now on we will restrict attention to fields which are invariant under simultaneous 
rotations and reflections in space and target space. Thus we assume that (fr = (4> a , </> 3 ) 
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with a = 1 , 2 is of the hedgehog form 

<t> a = sin g(r,t)k a , <f> 3 = cos g(r, t), (1) 

characterized by its topological charge k, defining the unit vector k a = (cos k9, sin k6) 
in ([I]) in terms of the azimuthal angle 9; and by the real (profile) function g of the 
polar coordinates and t which satisfies certain boundary conditions. The corresponding 
potential energy of the field (|l|) is 

1 f°° k 2 
E p = - (rg' 2 + -sin 2 g)dr, (2) 
4 Jo r 

where g' = dg/dr. This is normalized so that a static configuration has k energy. The 
boundary conditions are g(0,t) = n, in order to ensure a unique definition of 4> at the 
origin and g(r,t) ^ as r — > oo, so that E p converges. 

The standard Bogomol'nyi argument || is 

1 r°o, _ , fc 



< - / (Vrg r -\ psing) 2 (ir 

4 Jo vr 

/"OO 

= E p~2J ®r(cosg)dr 



= E p -k. (3) 

So the energy i£ p is bounded below by k; and E p equals k if and only if g' = —ksing/r, 
the solution of which is the static /c-soliton configuration 

g(r) = 2 arctan , (4) 

located at the origin, with a being a positive real constant which determines the soliton 
size. If a is large, the soliton configuration is flat and broad; while if a is small, it is tall 
but narrow. In fact, the height of the configuration (maximum of the energy density) is 
proportional to a~ 2 ^; while its radius (width) is proportional to a 1 ^ . Notice that, for 
k = the field is constant and the energy density is zero everywhere; while for k — 1 the 
configuration looks like a lump peaked at the origin; and for k > 1 it is a ring centered 
at the origin. In what follows, we will assume that in all cases k > 0, since taking k = 
does not test the ability of the model to handle nontrivial topologies. 

So far all we have done is to re-express the fc-soliton solution in terms of a real field 
g, which is a function of the polar radius r. It will now been seen how this description is 
useful in constructing discrete analogues of the Bogomol'nyi equations. 



From now on, r becomes a discrete variable, with lattice spacing h. So the real-valued 
field g(r, t) depends on the continuous variable t, and the discrete variable r = nh (n G Z, 
n > 0). The subscript + denotes forward shift, i.e. g + (r,t) = g(r + h,t) = g((n + l)h,t); 
and so the forward difference is given by Ag = (g + — g)/h. The question of how best 
to incorporate topological ideas into a lattice formulation has been the subject of much 
discussion, especially in lattice gauge theories contexts (cf. |§). One approach, following 
Speight and Ward is to begin with the same function cos g as appears in (|) and 



reconstruct the inequality, i.e. 

k(A cos g) = -D n F n , (5) 

where D n — > ^frg' and F n — > k sin g/y/r in the continuum limit h — > 0. The formula 
k(Acosg) = — 2k/hsm(2±p L ) sin(^ i ^) suggests the choices 

2f(h)^/n . (g+-g 



Vh V 2 

k ( g + + g\ 

: sin I — - — I , n > 0, (6) 



f(h)Vtoi V 2 

where f(h) is an arbitrary function of the lattice spacing with constraints f(h) ^ 1 as 



h — > and f(h) > yk/2 (see below). The origin must be treated in a special way since 
(|) are undefined when n = 0. One possibility is to arrange that D + F = identically. 
So choose 

[2k f g(h,t) \ 
D = -F = y x cos^^-j. (7) 

The potential energy of the lattice 0(3) sigma model field is defined to be 

h 00 
% n=0 

= fccos nsm ^—j+^mn 

As in the continuum case, it follows that E p is bounded below by k; and the minimum is 
attained if and only if D n + F n = 0. For models with different lattice spacing the effect 
of f{h) is to decrease the importance of the sin 2 g term in the energy density, although 
the total energy is still the same as in the continuum (k in our units). 

The kinetic energy can be defined by the simple choice 

E k = ^-^9\ (9) 

4 u=l 



where g = dg/dt. The boundary condition on g is that it should tend to zero at spatial 
infinity; this guarantees finite energy. For such fields, the total energy E t = E p + Ek 
is bounded below by k ; and this lower bound is attained if and only if g = 0, and 
D n + F n = for n > 0. [Recall that D + F = identically] 

This latter condition, i.e. D n + F n = 0, is called the Bogomol'nyi equation. It is 
a first-order difference equation, whose solutions (for the aforementioned boundary con- 
ditions) minimize the potential energy, and therefore, are also static solutions of the 
Euler-Lagrange equations 

J>iKnM) = -^, (10) 

since dE p /dg = at a minimum. So using the discrete Bogomol'nyi equations one gets 
first-order equations whose solutions are also static solutions of the second-order equations 
of motion. Moreover, these solutions have energy which is at its topological minimum 
value. 

The Bogomol'nyi equation D n + F n = 0, may also be written as 

g + 2f 2 n — k g 

tan T = 27^n tan 2' n>0 ' 



from which one sees that the function f(h) should be greater than yk/2 in order the 
profile function to be monotonic. So for k = 1,2 the form f(h) = 1 + h will work for all 



h; while for k > 3 we need a condition on h, i.e. h > Jk/2 — 1. 



The solution of (|TT| ) can be written down explicitly, i.e. 



( 7r, n = 0, 

g(nh) = \ (12) 
I 2 arctan^ii?™), n > 0, 

where 

r(n - k/2f) r(i + fc/2/ 2 ) 
n i> + A;/2/ 2 )r(l- A;/2/ 2 )' 1 J 

and z\ is an arbitrary positive constant which specifies, as in the continuum model, the soli- 
ton size. This is a static lattice fc-soliton solution located at the origin; which corresponds 
to a minimum of the energy in the k sector and thus, it is stable under perturbations 
which remain in that sector. From now on, we will concentrate on the k = 1, 2 sectors 
and therefore, we will take on the first one the function f(h) to be constant and equal to 
unity, for any h; while on the second one to be f(h) = 1 + h, for small h. 
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It would be nice to have a lattice analogue of the configuration width a 1//fc , which 
appeared in One possibility is to set 

a n = (nh) k tan (14) 

and then to define a = lim^oo a n , provided this limit exists. Indeed, a is proportional to 
Z\ h k . If E p = I]^L -E Pn , the energy density at the origin is E Po = k/(l + zf); therefore, 
E po is close to the Bogomol'nyi bound as z% — > 0. A diagram illustrating the profiles of 
the function g(nh) and the energy densities profiles are represented in figure 1 for k — 1 
and k = 2 with z\ = 15 and h = 0.19. 

The situation we wish to study is that of an isolated perturbed static fc-soliton con- 
figuration and investigate the effects of the perturbation. As it costs them no energy to 
shrink or expand they can shrink to almost a zero width (radius) configuration in the en- 
ergy density plot. Since the soliton configuration is described by a few points on a lattice 
it is difficult to decide what is meant by its width and how to calculate it. The lattice 
analogue will be a field configuration with g(0,t) = tt (due to the boundary conditions) 
and g(nh,t) = 0, for n > 0. In fact, this corresponds to a spike soliton (of almost zero 
width) in the continuum and thus, our interest lies in the study of the time dependence 
of the shrinking of a fc-soliton configuration. 

Since, there is no explicit solution in this case, one has to resort to approximation, or 
to numerical solutions of the equations of motion fllPf) , namely 

n _ 



1 k 2 
g = —[ksmg(h,t) + fsm(g(2h,t)-g(h,t))]-jj^sm(g(2h,t)+g(h,t)), 



ng = r [n ^_ 9) _ (B _ 1) ^_^ )] _^ [ ^ + - 5 ^ f ) ]i n>1 . 

(15) 



3. The Slow-Motion Approximation. 

There is a fundamental difference between the cases k = 1 and k > 1, which becomes 
apparent when one considers the slow-motion approximation, originally proposed in con- 
nection with monopole scattering |iTf . In this scheme one assumes that the field g is a 



static solution, but slightly perturbed. More precisely, since the energy is conserved, and 
due to the existence of the Bogomol'nyi bound, we may assume that a A;-soliton dynamics 



is obtained by restricting g to have the form of (12), with z\ now becoming a dynami- 
cal variable Z\(t). So the number of degrees of freedom is reduced from infinite to one. 
These static solutions form a manifold, which is equipped with a natural metric coming 
from the kinetic energy, and the evolution is given by the resulting geodesies. Since every 



configuration of the form (12) has the same potential energy, the kinetic energy may be 
taken as the Lagrangian; thus, the corresponding Euler-Lagrange equations are precisely 
the geodesic equations associated with the aforementioned metric. This approximation is 
a good one if the speeds are small (if Ej~ is small compared to E p = k). 

For k — 1, the requirement of finite kinetic energy means that z x should be independent 
of t at spatial infinity, so ruling out the slow-motion approximation. In other words, taking 
z\ to be a function only of t leads to a divergent kinetic energy. But when k > 1 there 



are sufficient powers of n in the denominator of (|12"D to keep the energy finite. For this 
case the slow-motion approximation has been considered in || [12j in order to study the 
dynamics of CP 1 lumps. However, they have not looked at the speed of the shrinking in 
any detailed. Let us concentrate on the k = 2 topological sector, where two solitons are 
sitting on top of each other at the origin, forming a ring structure. The Lagrangian is 

L = E k — E p 

= ^i)ia-2, (16) 

where 

^ ft \t<T^§r < 17 > 

The Euler-Lagrange equation of the system is 

2l(z 1 )z 1 + l'(z 1 )z 2 1 =0 - (l( Zl ) zf) = 0. (18) 



which may be reduced to quadratures: 



vt 



Zl (t) 



l{z x ) 



dz\ 



= AftCsO, (19) 

where Zi(0) = c, ii(0) = v. Recall that, Z \ determines the configuration size which evolves 
with t; more precisely, yfch is the initial width of the configuration and, v is the initial 
rate of change of the configuration width in each lattice site per unit time. In fact, v < 
corresponds to an initial contraction and v > to an initial expansion. 



The function A^Zi) decreasing or increasing depending on the value of Zi(t), which 
corresponds to contraction or expansion of the configuration. It is easily inverted to give 
the time variation of the configuration size, i.e. 

Zl {t)=K- h \vt). (20) 

Recall that, the lattice analogue of the configuration width is proportional to J~Z\h. In 
fact, the time taken for the configuration to shrink from the initial width to zero, i.e. to 
become a spike, is 

t c = -^-. 21 

v 

We believe that this analytical approximation is accurate for small \v\. 

The accuracy of the approximation has been tested numerically using a fully-explicit 
fourth-order Runge-Kutta algorithm with fixed time step 0.0053. The initial condition 
was a static 2-soliton profile whose width we perturbed to shrink with initial velocity 
0.1 lattice site per unit time (v = —O.lh). Simulations of duration 2985 time units were 
performed for h = 0.01. Inspection of the rate of change of z\(t) = tan(g(/i, t)/2) reveals 
close agreement with Z\{f) calculated from (p0|) (see figure 2). 

4. Dynamics of the Lattice 0(3) Solitons. 

Since the slow-motion approximation is expected to fail at high velocities (except for 
small h), we incorporate the notion of lattice solitons in a full numerical evolution scheme. 
Throughout the simulations, the extensive use of the difference equations ( [lSf) have not 
revealed any instabilities (the total energy is conserved). 

The lattice formulation necessarily has a spatial boundary at n — n max , say. Hence, 
the quantities that we are going to use in order to study the soliton dynamics will be 
calculated within some radius (n max ). Moreover, the infinite sums on the energies will be 
truncated. In fact, for k — 1 the finiteness of the grid imposes an artificial cutoff which 
provides a finiteness in the energy. On the boundary though, the fields are taken to be 
fixed in time, i.e. 

g((n max + 1) h, t) = g((n max + 1) h, 0), (22) 

since if one attempts to apply boundary conditions which allow the field to change with 
time at arbitrary distances, then the total energy of the system for k = 1, grows rapidly 
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and without bound. One may also, choose absorbing boundary conditions or may place 
the boundary far enough from the configuration (i.e., no radiation effects). But, in this 
scheme, the choice of the boundary conditions has no impact on the rate of shrinking. 

Moving on to the question of initial data, there are clearly many different types of 
perturbation which we could apply to the configuration, the only restriction being that 
we do not perturb the field close to the boundary. Since the evolution equations ( p~5D are 
second order the initial data must specify the field values g(nh, t) and its time derivatives 
g(nh,t) at t — 0. So the field configuration at t = is taken to be the static lattice one 
(|l2l) but slightly perturbed, i.e. 

2vZ 

g(nh,t)\ t=Q = 2 " 2 - (23) 

Physically the picture is this: there is a continuous interpolation between the inner region 
where v is the amplitude of the perturbation (as in the slow-motion approximation) and 
the outer one where there is no perturbation at all. This class of perturbation reveals all 
the qualitative types of behaviour that can occur. 

So we have a fc-soliton configuration whose centre remains fixed, but whose radius 
decreases to a minimum (close to zero) and then increases again. More precisely, the 
initial perturbation (for v < 0) tends to shrink the configuration, while large burst of 
radiation travel outwards at the speed of light (see figure 3), together with a residual 
motion in the central region occupied by the soliton. When the radiation reaches the 
boundary, is reflected back, reabsorbed by the configuration which expands and then 
another pulse is emitted a short time later; and so the process repeats. The data were 
produced by the aforementioned Runge-Kutta algorithm for v = —0.1 and z\ — 1, on a 
lattice of unit spacing (h = 1) in the k — 1 sector. 

We are interested in the speed of the shrinking in detail. In order to analyze the results 
of the numerical simulations we look at the dynamical quantity g(h,t). Since we are on 
a lattice, the soliton will be highly localized (a spike) when the profile function at the 
first site (n = 1) and consequently, at all others (n > 1) will be zero. Then, the soliton 
configuration occupies essentially only one lattice site, while the slope of the curve g(h, t) 
correspond to the power law of shrinking. 

In figure 4 we present the time dependence of the field g(h, t) for a single soliton and a 
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2-soliton configuration. The results are derived from a relatively small mesh n max = 200 
(in fact, they do not change for larger mesh sizes), for v = —0.1 and z\ = 1. In the single 
soliton case, we make the simple choice f(h) = 1 whereas h — 1; while in the 2-soliton 
case, we take f(h) — 1 + h with h = 0.01. (Note that, figure 4(b) corresponds to figure 
2.) The field configuration saturates the Bogomol'nyi bound throughout the numerical 
evolution and due to that in the k — 1 sector the lattice spacing is comparable to the 
soliton size without compromising its behaviour. As it is clear for the graphs the curves 
are nearly straight, confirming the power law for the rate of shrinking. The linear curve 
in figure 4(b) is due to the fact that the lattice spacing is small compare to the size of the 
topological soliton and thus, the model is closer to the continuum one. Let us conclude 
with the observation that our results are consistent with the ones obtained by studying 
the continuum 0(3) sigma model (cf. ||). 

5. Conclusions. 

We have studied the time evolution of the lattice 0(3) sigma solitons which were 
allowed to shrink. We have used two different methods to analyze numerically the soliton 
shrinking. The first method is based on the slow-motion approximation and leads to 
an ordinary differential equation. The second one consists of integrating numerically 
the (l+l)-dimensional semi-discrete equations. From all the numerical studies we have 
performed we conclude that a single soliton and a 2-soliton configuration of the 0(3) 
lattice sigma model does shrink to a highly localized (spike) lattice soliton, linearly with 
time. If one is close to the continuum limit, in a sense that the lattice spacing is small 
compared to the size of the topological solitons, then there may not be much difference 
between various lattice versions of the continuum system. An advantage of the lattice 
model described in this paper is that the lattice spacing can be relatively large without 
compromising the soliton dynamics. 

Let us conclude by stressing once again that the slow-motion approximation works 
very well. 
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Figure 1: (a) Profiles of g for charges k = 1,2. (b) Profiles of the energy densities E Pn /(2imh) 
for the charges in (a). The k = 2 energy density is ring shaped. 
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Figure 3: Radiation emitted by the 1-soliton solution. 
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Figure 4: The variation of g(h, t) over the range (a) < t < 12.6 for a slowly shrinking 1-soliton 
lump and (b) < t < 16 for a slowly shrinking 2-soliton ring, for the numerical evolution. 
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